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ABSTRACT 

We analyse the KSB method to estimate gravitational shear from surface-brightness mo- 
ments of small and noisy galaxy images. We identify three potentially problematic assump- 
tions. These are: (1) While gravitational shear must be estimated from averaged galaxy im- 
ages, KSB derives a shear estimate from each individual image and then takes the average. 
Since the two operations do not commute, KSB gives biased results. (2) KSB implicitly as- 
sumes that galaxy ellipticities are small, while weak gravitational lensing assures only that 
the change in ellipticity due to the shear is small. (3) KSB does not invert the convolution 
with the point-spread function, but gives an approximate PSF correction which - even for 
a circular PSF - holds only in the limit of circular sources. The effects of assumptions (2) 
and (3) partially counter-act in a way dependent on the width of the weight function and of 
the PSF. We quantitatively demonstrate the biases due to all assumptions, extend the KSB 
approach consistently to third order in the shear and ellipticity and show that this extension 
lowers the biases substantially. The issue of proper PSF deconvolution will be addressed in a 
forthcoming paper. 

Key words: gravitational lensing methods: data analysis cosmology: observations 



1 INTRODUCTION 

Cosmic shear measurements are a potentially powerful probe of 
structure growth at intermediate to late epochs of the cosmic his- 
tory. Sufficiently precise measurements should be able to con- 
strain both the amount, distribution and fluctuation amplitude of 
dark matter and the time evolution of dark energy ("see lSartelmannl 
l201Ct for a recent review). Signi ficant cosmic -shear signals have 
been detected in many studies (see lSacon et aU2 000: Kaiser 200ol : 
ISenjamin et al ] |2007L for recent examples) and evidence for ac- 
celerated expansion Ischrabback et al.l(2009t) has also been found. 



Cosmic-shear measurements are the primary motivation for sev- 
eral dedicated surveys that are proposed or upcoming (e.g. EU- 
CLIeQ, JDEIV0, DE^l LSStQ). They will greatly increase the sur- 
vey area and the number of observed galaxies, and hence lower the 
statistical uncertainty of the shear estimates. Howev er, the anal- 
ysis of synthetic data in the GREAT08 challenge jBridle et al.l 
12009) shows that the accuracy of shear estimation methods is cur- 
rently not sufficient to fu lly exploit the next-generation surveys 
jAmara & Refregiej2008h . 

Particularly concerning are systematic biases in shear esti- 
mates, which do not vanish when averaged over a large ensem- 
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ble of lensed galaxies. These biases often stem from assump- 
tions made in the derivation or implementation of shear estimation 
methods, which do not hold in reality. For instance, the models 
used to describ e the galactic shape ma y systematicall y differ from 
the tr ue shape jLewi j|2009l : IVoigt & Bridle 2010; M elchior etakl 
I2OIOI) . Model-independent approaches may therefore be favored as 
they are not limited by the peculiar ities of an underly ing model. 
The widely employed KSB method jKaiser et alj|l995h is model- 
independent because it expresses the lensing-induced shape change 
by combinations of moments of the galactic light distribution. 
However, it relies on several assumptions regarding (1) the strength 
of the apparent distortion, (2) the width of the window function, (3) 
the mapping between convolved and unconvolved ellipticity, and 
(4) the ellipticity of the PSF. 

It has been noted that the accuracy of the K SB method has 
a problematic dependence on items (2) and (4) jKuiikenlll999l : 
lErben et alj2001 t lHoekstra et al.ll998t) . In this work we investigate 
the reasons for problems encountered with shear estimates from 
KSB by a rederivation of its fundamental relations. We take par- 
ticular care of inspecting the assumptions made, and show if and 
how improvements to the original KSB relations can be incorpo- 
rated such that the shear estimates remain free of bias in a wider 
range of galactic and PSF parameters. 

In Sect. 2, we briefly review the basic relations of gravitational 
lensing and show how to construct a shear estimate from observed 
image ellipticities. In Sect. 3, we introduce and test three variants 
of KSB based on a linearized relation between shear and ellipticity. 
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and a novel one which employs a third-order relation. We inves- 
tigate the validity of the PSF-correction approach in Sect. 4 and 
comment on the possibility of an improved correction for PSF el- 
lipticities. We conclude in Sect. 5. 



2 WEAK LENSING BASICS 

This section summarises the basic weak-lensing concepts that 
will be used later. Fo r a c omplete overview we refer to 
iBartelmann & Schneided j20011) . An isolated lens with surface 
mass density has the lensing potential 

4G DiD, 



Di: 



■J (f9'E{e')ln\9-9'\, (1) 



where G and c are the usual constants and D^s.is are the angular- 
diameter distances between the observer and the lens, the observer 
and the source, and the lens and the source, respectively. 

To sufficient accuracy, light rays are deflected by the angle 



a(6») = V*(6») 



(2) 



which relates the angular positions of the source /3 and the image 
on the sky by the lens equation 



(3 = -a{0) 



(3) 



If the lens mapping changes little across the solid angle of a source, 
the lens mapping can be locally linearised to describe the image 
distortion the Jacobian matrix 

4 = _ / _ \ _ / 1 - ^ - 71 
^- dQ- y"'' de.dOj )-\ -72 
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with the convergence 



( ISchneider&Seitzlll995h . where is the unlensed (intrinsic) el- 
lipticity. This relation holds as long as the lens mapping can be 
locally linearised. Information on the intrinsic ellipticity of a single 
object is not accessible. Reasonable shear estimates thus require av- 
eraging over many galaxies in a region where g can be considered 
constant, assuming that the average of vanishes. 
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If the coordinate frame is rotated such that only one shear compo- 
nent does not vanish, Eq. I ll2t is solved by 



(x) 



2(1 



+ 



ix)' 1 



) 8 (l-al) 



'^ + 0{{xf) 



(13) 



where is the standard deviation of the intrinsic ellipticity distri- 
bution. In the derivation of the equation above we neglected higher 
order moments of the intrinsic ellipticity distribution. Note that the 
average ellipticity appears in this equation, and that the relation be- 
tween the average ellipticity and the shear is generally non-linear. 
We recall here that other ellipticity estimators can be defined in ad- 
dition to the one presented in Eq. dlOl l. Another common estimator 
is 



(Qll-Q22)+2iQl2 
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w hich has a perfect respon se to shear, i.e. the shear responsivity is 
1 l lSeitz&Schneideilll997l) . However this estimator is considered 
more noisy and therefore not commonly used in weak-lensing mea- 
surements, and in particular it is not used by KSB. For this reason 
we will employ x ellipticity estimator rather than e throughout 
this work. 



and the two components 
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of the complex shear 7 = 71 + 172. Image distortions measure the 
reduced shear 



g = 



7 



(7) 



1 - K 

instead of the shear 7 itself. To linear order, G and (3 are related by 

ft = A.jO' . (8) 



2.1 Shear estimation 

The shape of an extended source can be descibed by angular mo- 
ments of its surface brightness distribution I{9), 



Qij . 
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Q is the total flux, Qi defines the centroid of the image, and higher- 
order moments provide information on the image's morphology. 
Combinations of second moments are used to quantify the image's 
ellipticity, which we introduce as 

(Qll -Q22)+2iQl2 



'''' Qii + Q22 

The complex ellipticity x is related to the reduced shear g by 



X 



l + |gP-2J?(3x*) ' 
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3 SHEAR MEASUREMENTS 

In practice, shear estimates are obtained from small and noisy back- 
ground galaxies. The observed shape of any object is the result of a 
convolution of its intrinsic surface brightness with the point 

spread function P(0). The convolution tends to make the object 
more circular or to imprint a spurious ellipticity on it if the PSF is 
not isotropic. Moreover, any measurement of moments has to incor- 
porate a weight function in order to suppress the pixel noise dom- 
inating at large spatial scales. Convolution and weighting change 
the surface brightness to 



r'"'{0) = w{0) J i°(e')P(e - e')cFe'. 



(15) 



Since we are interested in the object's unconvolved and unweighted 
shape, we need to correct these two effects. In this Section, we first 
assume P{6 — 0') — >■ S{6 — 0'), i.e. we neglect the PSF convolu- 
tion, and postpone the PSF correction to the following Section. 



3.1 Standard KSB 

We review in this section the standard KSB formalism, neglecting 
PSF convolution. In this situation the only complication is given 
by the presence of the weighting function for the computation of 
moments, which modifies the relation between shear and ellipticity 
given by Eq. dl lb . Weighting changes Eq. dlOb to 



Tr{Q) 



(16) 



KSB: biases and corrections 3 



with 



ria 



el - el if a = 1 

26I16I2 if a = 2 



(17) 



Note that also TrQ in Eq. l ll6t is evaluated using weighted mo- 
ments. Using Eq. ((Sj and the conservation of the surface bright- 
ness, I°'"'{0) = P{A0), we can infer the surface brightness in the 
source plane. From its second moments. 



= {det A) A.kAa j d'er'"' {6)6^61 
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we form the ellipticity 
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where 



TriQ^) 



(19) 



(20) 



^^^no.-2g^\e\-' + {-irn^{gl-gl) + 2g^g2ni ■ (21) 

The relation between the two filter scales in Eq. dlSt is given by 
o-^ = (l-fi:)^(l + |5|^)cr^ andthemultiplicativeterm (dety4)(l- 
k)^ in Eq. ( I20t will cancel out once Tr{Q") is written in terms 
of Tr{Q). Note that Einstein's sum convention is not implied in 
{-l)°'rja, and that 



t _ r 7,2 if a = 1 
\ ryi if Q = 2 

We adopt this notation for a general tensor. 



fit 



fi2/3...C if Q = 1 

fii/3...C if a = 2 



(22) 
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Combining Eqs. l ll6b and l |19t gives a more complicated rela- 
tion between ellipticities in the source and in the lens planes than 
Eq. Jl lb due to the presence of the weight function. Keeping only 
first-order terms in g, this relation is 



Xa — Xa — 9 t^ap^ 



jKaiser et alJl995l : lHoekstra et al.ll 19981) with 



p:.p = -2 



'Tr{Q) 
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2S^f, (25) 



and 
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The notation we use here follows lBartelmann & Schneideil ( 1200 ll) . 



3.2 Shear estimates 

Equation i24l directly relates the measured weighted ellipticity x 
to the shear g if the intrinsic ellipticity of the source x" is known. 



Since x and x" cannot be disentangled for individual galaxies, av- 
erages over ensembles of images are necessary to estimate g, 

iXc) - iXc) = (g^P^p) (5a) = Pa = (P-'^Vpix^) . (27) 
However, the quantity which is commonly computed is 

{9c.) = {{P-'-Ux^) , (28) 



lErben et alj ( 1200 ih , assuming that ((P°'')"^x''> = 0. This condi- 
tion is not guaranteed since P"'* itself depends on x- The symbol 
ga denotes the shear estimate obtained by solving Eq. l l24t with 
X" = 0. We introduce it since g is not the true shear (which is inac- 
cesible for a single galaxy), but the shear one would measure if the 
source was circular. The true shear g is then sought by averaging g. 
Equation dl It shows that for x" = and W{x) = 1, (g) is related 
toxby 



(5> = 



l^^/^^ 

X 
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(29) 



In general, (g) differs from the true shear g computed in Eq. l ll3t . 
Assuming that g <C 1, meaning (x) ^ 1, and the distribution of 
the intrinsic ellipticities to be Gaussian with standard deviation a^, 
the difference can be written as: 



from which 



9-r9}^^4 
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For a realistic ~ 0.3, the bias introduced by averaging shear 
estimates instead of ellipticities is ~ 2%. 

Moreover, averaging shear estimates does not allow one to as- 
sume that g is small, as done in the original derivation of P^'\ since 
it is always of the same magnitude as x- In coordinates rotated such 
that g has only one non- vanishing component, and in absence of a 
weight function, the relation between x and g provided by KSB in 
Eq. ^ is 



-KSB ^, Xl , Xl , Xl , 

^ ~2 2 2 



(32) 



Obviously, this is correct only to lowest order. Comparing Eqs. 
i29\ and l |32t , the error made by KSB in the shear estimation 
is a function of the measured ellipticity and scales as (3x^/8 + 
7x^/16). Typically, |x| G [0.5. ..0.8], implying that the bias KSB 
introduces in the shear estimate (without weight function) is in 
the range [6... 33]%. The reason for this bias comes from the fact 
that second- or higher-order terms in g have been neglected in the 
derivation of Eq. ( 124b . while terms like x^g have been kept. Once 
g is identified with g, these mixed terms are effectively of the same 
order as the terms. In a consistent first-order relation between x 
and g, only the first-order term in xg can be considered. Then, P^'^ 
looks like 

2Ba 



psh,(0) 



3a/3 

TrQ 



+ 2&aP 



(33) 



We shall refer to this approximation as KSB 1 . In this case, the so- 
lution for g is 



g 



X 
2 



(34) 



The error on the shear estimate made by KSBl scales like 
(— X'VS — X^/16), leading to an underestimate which is consid- 
erably smaller than the overestimate given by KSB. However, as 
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discussed before, x is practically never small, meaning that first- 
order approximations may be poor. 

In a frequently used variant of KSB, is approximated by 
half its trace (KSBtr hereafter). 



(35) 



This is usually justified saying that the trace is less noisy than the 
inverse of the full tensor, as we shall show in Sect. 3.3. This state- 
ment is certainly correct for large ellipticities. However, it turns 
out to work much better than the full tensor even in the absence of 
noise, PSF and weighting. The reason is that it leads to the relation 



3 5 
-KSBtr ^, X . X . X , 

^ ~2 4 8 



(36) 



between g and x, which biases the shear estimate by x^/8+x^/16. 
We can summarise the preceding discussion as follows: 

(i) KSB incorrectly approximates Eq. i ll It : 

(ii) KSB 1 is mathematically consistent; 

(iii) KSBtr approximates Eq. dllb better even though it lacks 
mathematical justification; 

(iv) No KSB variant discussed so far is correct to third order in 

X- 



3.3 Third-order relation between g and x 

We now derive a consistent third-order relation between x ™d g, 
including the effects of the weight function. We fol l ow cl osely 
the approach in Sect. 4.6.2 of lBartelmann & SchneideJ jlOOlh , and 
use Einstein's sum convention. We start from Eq. J19t and Taylor- 
expand the weight function around g = Qto third order in g. 
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(37) 



(38) 



was used. Note that the derivatives of the weight function are taken 
with respect to 6^ . Truncating the series at a given order implies 
that the final result will depend on the shape of the weight function. 
We proceed with the calculation of XaTr{Q'') to third order 



X^TriQ^ 



(detA)(l - k)2 



XcTriQ) - 2g^Bc,p + 2g^g^Dc,p^ ~ 2g^Tr{Q) + 'ig^g^Lp + 

- 4Kp.,g^g^g^ + {-lY'igl ~ g!)xcTr{Q) 

- 2{~irB^^{gl - gl)g^ + 2gig2xiTr{Q) + 

- ^Bl^g^g^g^ - \u^Mi9^9''9' + 

(39) 



where the definitions 

Dq/37 = ^ / (f 91°^^ {e)W" r^aTipri^ , 



(40) 



appear. La and B^p are given in Eq.[26] In the same way, we eval- 
uate 

(det^)(l - 



Tr{Q){l + \g\') - 2g"L, + 2p"/7^,^ 



'^g°'XaTr{Q) + Ag°'g^Bai3 - 'IDap^g'^ g'^ g^' - -^Jai3',g°'g^ g^ 
= Tr{Q){l + f{g)), 

(41) 

where we implicity defined f{g) and 

= ^ f d^oi°'"'{0)w"'\ef-navpn-y ■ (42) 

From these quantities, we compute 

^ ^ XcTr{Q){l + f{g))-xiTr{Qn ,,,, 

^■'^ Tr{Q){l + f{g)) ' 

This equation holds exactly in absence of a weight function. If a 
weight function is included, f(g) is at most of order 0.02, and we 
shall consider Eq. ( 143 1 exact to third order. After some algebra we 
find 

/[Pa/3 + g^iRafi-, + g^Scp^s)] + H„ + 



Xc - Xc. 
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where 
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and 



Xcc 



TrQ 



2gw. [x^ - j 



,2 / ,-.a, 2 2s I 2Bafig' 
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(46) 



Introducing second and third-order terms leads to a non-linear rela- 
tion between x ™d g which needs to be solved numerically. More- 
over, sixth-order moments of the observed surface-brightness dis- 
tribution appear in Rap-f and eigth-order moments in Uap-iS be- 
cause of the Taylor expansion of the weight function to third order. 
We discuss in the following Section how to deal with the non-linear 
relation between shear and ellipticity and possible noise issues due 
to the appearance of higher moments. 

3.4 Tests 

We now show the results of simple tests carried out to check how 
well the four variants of KSB estimate the shear. We consider a 
circular source (x" = 0) with a Sersic brightness profile. 



7(r) = Jo exp 



(47) 
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Figure 1. Shear estimate gi as a function of the applied shear for noise-free and unconvolved Sersic-type galaxy images as provided by KSB (red line), 
KSBtr (green line), KSB3 (blue line), and KSBl (magenta line). In the left panel no weighting function has been used to measure moments of the light 
distribution, while in the right panel a Gaussian weighting function has been employed with a width equal to the size of the object. The effective galaxy 
radius was iJe = 2 pixel, the Sersic model was tenfold oversampled, and the image sidelength was 40 pixels. 



where is the radius containing half of the flux and the Sersic 
index and is a constant which depends on ris. This type of pro- 
file is identical to a Gaussian for Us — 0.5 and is steeper in the cen- 
tre for Us > 0.5. In the following test, we assume rij, — 1.5, which 
represents t he averag e value for rather bright galaxies in the COS- 
MOS field jSargent et al. 2007). We shear this profile by a variable 
amount gi, keeping g2 — 0.1 fixed, using Eq. ([Sj. For all follow- 
ing tests, the effective galaxy radius was Re = 2 pixels, the Sersic 
model was tenfold oversampled, and the image sidelength was set 
such as to not truncate the galaxy at the image boundary. Then, 
we measure the ellipticity as defined in Eq. JlOb . Since the model 
galaxy is intrinsically circular, the source ellipticity is entirely gen- 
erated by the applied shear which is varied in a wide range such 
as to mimic the intrinsic ellipticity dispersion. The weight function 
has been chosen as Gaussian with a = 2i?e. We repeated this test 
assuming a flat weight function (W(x) — 1) in order to estimate 
how much the different approximations in deriving P"'^ affect the 
measurement. The results are shown in Fig.[T] 

In absence of a weight function (left panel of Fig.O, the per- 
formance of the four variants closely follows the analytic behaviour 
worked out in Sects. 3.1 and 3.2: KSB severely overestimates the 
shear for large gi, while KSBtr and KSBl better approximate the 
shear. KSB3 returns the correct shear under this condition. 

The weight function renders the image more circular and thus 
reduces the measured x- This means that the high-order terms in x 
contribute less to the shear estimate. Therefore, the deviation from 
the correct result is significantly lower for all the methods (right 
panel of Fig.lTJ. This is not true for KSBl, which allows only a 
first-order correction for the weight function. 

We also investigate the behaviour of the four KSB variants for 
realistic pixel noise. The average result for 100 galaxies is shown 
in Fig. [2] 

KSBtr is the only method for which no matrix inversion is re- 
quired. It is thus not surprising that it exhibits the lowest standard 
deviation for all values of g. KSB and KSB3 have a comparable 
amount of noise even though KSB3 involves the computation of 6th 
and 8th moments of the light distribution. The reason is that these 
higher-order moments are computed using the second and third 
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9l, true 

Figure 2. Shear estimate gi as a function of the applied shear for noisy 
but unconvolved Sersic-type galaxy images as provided by KSB (red line), 
KSBtr (green line), and KSB3 (blue line). The total flux of the source was 
fixed to unity, the noise rms to 10""^. The average is taken over 100 objects. 
Errorbars denote standard deviation of the mean. 

derivatives of the weight function. There is no price to be payed (in 
terms of measurement noise) in using KSB3 instead of the simple 
KSB description. We also investigate how much the measurement 
of one component of the shear is affected by the value of the other 
component. For this case, we also studied the case of unweighted 
and weighted moment measurements. The result is shown in Fig. [3] 
The obvious cross-talk between the two components is not surpris- 
ing for KSB, KSBtr or KSBl since all terms which mix gi and g2 
were neglected in the calculation. Introducing third-order correc- 
tions, the estimate of one shear component becomes almost inde- 
pendent of the other component. 

Finally, we study how much the bias in the shear measure- 
ment depends on the width a of the weight function W. We vary 
the width within [2Re, oo). The result is shown in Fig. |4] KSB 
and KSB 1 exhibit a strong dependence on a, while KSBtr is more 
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Figure 3. Shear estimation cross-talk for 32 as a function of tlie applied shear g\ for noise-free and unconvolved Sersic-type galaxy images as provided 
by KSB (red line), KSBtr (green line), KSB3 (blue line), and KSBl (magenta line). In the left panel no weighting function has been used to measure 
moments of the light distribution, while in the right panel a Gaussian weighting function has been employed with a width equal to the size of the object. 




Figure 4. Dependence of the shear estimate gi on the size of the weight- 
ing function width as a function of the applied shear for noise-free and 
unconvolved Sersic-type galaxy images as provided by KSB (red), KS- 
Btr (green), KSB3 (blue), and KSBl (magenta). We consider a Gaussian 
weighting function with width a = [2/?^, 00]. The lower limits corre- 
spond to (T = 2Re, and the upper limits are identical to the unweighted 
case shown in Fig.[T] 



robust, and KSB3 is almost independent of cr. Due to the poor cor- 
rection of the weight-induced change of x, KSBl performs poorest 
in this test. For KSB, the reduction of % due to the weighting limits 
its strong non-linear response such that the bias decreases for nar- 
row weight functions. As KSB3 employs the best description of the 
weighting, it performs excellently in this test. 

In all tests carried out so far, we have assumed that the in- 
trinsic ellipticity of the object vanishes, x° = 0. This is of 
course idealised since galaxies have an intrinsic ellipticity disper- 
sion. In order to test the performance of the four methods for an 
isotr opic source-eUipticity d istribution, we apply the so-called ring 
test jNakaiima & Bemstein|[2007i") . We construct an ensemble of 
test galaxies falling on a circle in the ellipticity plane, shear them, 
measure their shapes, and take the mean. We choose an intrinsic 
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Figure 5. Shear estimates of a sample of Sersic-type galaxies with x" = 
0.3 after shear g = (0.1, 0.05) is applied. Red dots are the results from 
KSB, magenta dots from KSBl, blue dots from KSBtr, and green dots from 
KSB3. The dots in the center show the position of the ensemble averages 
of the estimates. A zoom of the central region is shown in the small panel, 
where the intersection of the dotted hnes indicates the outcome of a perfect 
measurement. 



ellipticity \x''\ = 0.3 and apply the shear g = (0.1,0.05). The 
result is shown in Fig. [5] A perfect method would recover the cor- 
rect shear after averaging over all test galaxies. Not surprisingly, 
we find that KSB is unable to recover the correct shear from the 
averaged individual shear estimates since they depend non-linearly 
on X- This leads to an average overestimate of « 35 % if the shear 
is aligned with the intrinsic ellipticity. As the other variants have 
lower non-linear error in the x-g relation, the mean values are bi- 
ased by ^ 20 % (KSBl), f» 5 % (KSBtr) and ^ 1 % (KSB3). 
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4 PSF CONVOLUTION 

Any measured galaxy's ellipticity is the result of three distinct 
physical processes: intrinsic ellipticity, lensing, and PSF convolu- 
tion. As discussed above, lensing maps the galaxy's light distri- 
bution from the source to the lens plane, distorting its shape. The 
relation between galaxy ellipticity and shear can be found solving 
Eq. dllt if there is no weight function, or Eq. l l44t if weighted mo- 
ments are used to define the ellipticity. In general, the equation one 
needs to solve to relate ellipticity to a shear estimator g has the 
implicit form 



X = fi9,x) 



(48) 



On the other hand, PSF deconvolution maps the observed ellipticity 
from the image plane (on which the object is lensed and convolved 
with the PSF) to the lens plane (on which the object is lensed only) 



(49) 



Thus, the relation between observed ellipticity and the shear esti 
mator in presence of PSF convolution is the solution of 

X""' ^h-'[f{g,h{x"'n)]■ 



(50) 



If the PSF is perfectly circular, the only effect of /i is a cir- 
cularisation of the object, otherwise the PSF induces additional 
anisotropic distortions. Therefore, it is crucial to properly correct 
these two effects in order to realiably estimate the shear. We con- 
sider in the following the case of a spherical PSF and briefly discuss 
the case of an anisotropic PSF in Sect. 14. 11 

Instead of carrying out a proper PSF deconvolution first and 
then estimating the shear using the unconvolved ellipticity, as sum- 
marised by Eq. l |50l l, KSB links the observed ellipticity to the shear 
by the following approach: 



obs sh g 

Xa Xa XoL ? 



(51) 



where x'^ is given by Eq. i44i . and x^ is 

xi=Pl7iP'"''l-'pXt'*- (52) 
P"™ is the so-called "smear polarisability tensor" and has the form 



a/3 



where 



Tr{Q) 



(53) ■ 



M = J dPei{e)w (^J|f ^ , 

= ^ J d'ei{0)W' (^1^^ 77a, and 

= ^1 <fei{e)w" {^^^ n.vp- 

L'a has to be interpreted as La calculated with the second deriva- 
tive of the weight function, while Tr{Q') and Tr{Q) are cal- 
culated with the first derivative of the weight function. We refer 



(54) 



to Sect. 4.6.2 of iBartelmann & Schneidej ( 1200 ih for a complete 
derivation of Eq.|51| 

Since x"™ encodes the action of lensing (cf. last Section on 
the appropriate forms of this mapping), we can rewrite Eq.|51|as 



X 



fia V ) - yHa v ) 



(55) 



It is important to note that the lensing-induced mapping is now 
evaluated with the observed, i.e. convolved, ellipticity instead of 
the unconvolved ellipticity. This approach therefore requires the 



correction term x^ , which corresponds to a correct treatment of 
the PSF convolution (Eq. l lSOb ) if and only if 



a i ~ obs \ 

X(9,X ) 



f(g,x°''l 



h-\f{g,h{x° 



(56) 



We study now a very simple but instructive case. We assume a 
perfectly circular source, no weight function, an isotropic PSF, and 
shear oriented in a single direction. Then, P^™ becomes diagonal. 



M 



(57) 



In a forthcoming paper (Melchior et al., in prep.) we shall demon- 
strate how to do a proper PSF deconvolution, using the moments of 
the PSF and the convolved object, and show that the mapping h be- 
tween the convolved ellipticity x"*"' and the unconvolved ellipticity 
X in the lens plane is given by 



Kx° 



i~ A{x°^''y 



(58) 



where 



A ■ 



M TrQ' 
TrQ M* 



(59) 



is a functionof the observed ellipticity (as shown in Fig. ^ and of 
the size of the PSF (as shown in Fig.|7]l, and is bound to [0, 1] . If 
the shear has a single component and there is no weight function 
involved in the measurement, /(gt, x"*"*) is 



J'yy^X ) l + g2_ 2g^obs 

According to Eq. i52l . in the KSB formalism has the form 



(60) 



x'{9,x'"'n = A{x''nfr9,o)- 



(61) 



,X°'n = 2~gA{x'"'n 



In particular, in standard KSB, x^{ 

Substituting this expression for x^ in the Ihs of Eqs. ( 156b . we 
can conclude that KSB gives a proper description of PSF decon- 
volution only if the function f{g, x°''") can be decomposed into a 
product of two functions, one depending on g only and one on x"'" 
only. This is by no means guaranteed. A detailed analysis reveals 
that there are two limiting cases in which Eqs. ( 1561 ) holds: 

• The PSF width vanishes: 

A(x°'n=x'ix'"'l = O^X'"' =X- 

• The observed ellipticity vanishes. 

While the first case is trivial (but irrelevant), the second case 
can only be realized - for any finite PSF width - by a conspiracy of 
intrinsic and lensing-induced ellipticity. 

To study in detail the error commited by KSB in the attempts 
to correct for the PSF convolution, we solve Eq. ( |51t explicitly, 
employing the four variants x"'^ of mapping x onto g presented in 
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the previous Section, 



-KSB 



f + f A'(0) + 



(62) 



Xo 



[(l-yl(0))(l + A"(0)/2)+A'(0)'] 



+ o(xl) 

-KSBl _ X 

*' ^ 2 

xo , xd 



KSBtr ^ — + — y4.'(0) + 

+ ^[(1 -yl(0))(l + A"(0)) + 2A'(0)'] 
+ O(xo) 

~3 

+ + 4^'(0)2 + 2A"(0) - 2A(0)(2 + A"(0))] 



+ 0(Xo) 



where 



Xo 



X 



1-^(0) 



(63) 



and A'(0) and A" {Q) are the first and the second derivatives of 
A(x ) computed for x°^" = 0- If the PSF correction works per- 
fectly, the relation between x and g has the same form as the exact 
unconvolved solution of Eq. l |13l l, 

-,3 



9 ■ 



1 + Y + o(x') 



(64) 



We note first of all that Eq. l l62b is written in terms of xo, while 
Eq. ( I64t is written in terms of x, meaning that in general the solu- 
tions are different already at first order. However the error at first 
order {x°*"' ^ 1) is mostly of order 10~^ and therefore negli- 
gible. In the limit of a very wide PSF A(x) — A(0) we find the 
deviations from the exact solution b = g — g'^^^ ". 



(65) 



b"^"" = i^x^ + o(x^) 

It is worth noting that the PSF correction introduces a bias with 
preferred direction: Shear estimates decrease as the PSF width in- 
creases. 



4.1 PSF anisotropy 

An anisotropic PSF introduces spurious ellipticity in the image 
plane which must be corrected. The appropriate correction in KSB 
relies on the hypothesis that the PSF can be considered almost 
isotropic. This enables its decomposition into an isotropic part p"° 
and an anisotropic part q, 



P{6) 



(66) 



Even this dec omposition can be problematic for certain PSFs 
jKuiikenl 19991) . For example, a PSF given by the sum of two Gaus- 
sians with constant ellipticity does not fulfill the equation above. 



Assuming that Eq. l l66b is valid, one can find a relation, valid to 
first order in q, between the observed and the isotropic ellipticity. 



iso obs / j-tsm\ 8 

Xa =Xa - (Pcfi )q 



(67) 



The term qa , carrying information on anisotropies in the PSF, can 
be determined from the shape of stars using the fact that their 
isotropically smeared images have zero ellipticity (x*'"*" ~ 0)> 

Qc. = [P )c,f5Xfi (68) 

Once q has been determined, we can use Eq. l |67| ) to compute the 
isotropic from the observed elli pticity. For a detailed calculation 
we refer again to jBartelmann & Schneider 2001). In the deriva- 
tion, all the terms containing moments of q higher than the second 
have been neglected as well as quadratic and higher-order terms in 
qij. If one wants to extend this calculation to higher orders in q, 
derivatives of the observed surface brightness 1°^" would appear in 
the calculation because the assumed equality to V^" (hypothetical 
surface brightness for vanishing q) does not hold any more. The 
fact that derivatives of the observed profile need to be considered 
renders it practically impossible to incorporate higher-order correc- 
tions for q, since 1°'^" is a noisy quantity. This means that within 
the KSB framework it is not possible to correct properly for highly 
elliptical PSFs. 

If the determination of q is wrong, so is the estimate of x"° 
(the ellipticity of 7*"°), and the error will propagate to the final 
shear measurement in an almost unpredictable way jErben et"al] 
l200lHKuiikerilll999l) . This could happen if the anisotropy of the 
PSF is too large for a linear treatment, or if the PSF cannot be de- 
composed into an isotropic and an anisotropic part. 

4.2 Tests 

We perform the same tests as in the previous Section, but with an 
additional convolution with a Moffat- shaped PSF, 



P(r) = (1 + ar'^y 



where 



>l//3-l 



(69) 



(70) 



(FWHM/2)2 

controls the size of the PSF and /? regulates its steepness. In or- 
der to ensure vanishing flux at large radii, the PSF is truncated at 
5 FWHM, and the appropriate value at that position is subtracted 
from P{r). 

We begin studying the case of a flat weight function, W{x) = 
1, for which we derive the behaviour of the four KSB variants in 
Eq. l l62t . The key quantity for describing a spherical PSF is given 
by A as defined in Eq. l l59t . which is a function mainly of the PSF 
width and mildly of its steepness for a given galaxy (see Fig.|7]l. 
We investigate the performance of the four methods as a function of 
the shear for a fixed PSF width. We choose FWHM = 0.5 iie and 
/3 = 2 to mimic a space-based observation, and FWHM — 5 
and P — 5 to mimic a ground-based observation. The results are 
shown in Fig.|9] In the first case, KSB3 gives the best result, while 
KSBtr is the best approximation in the second case, as expected 
from Eq. l l65t . 

We next investigate the response of the four methods to the 
size of the PSF for a given g (Fig. [8). We choose gi — 0.4 and 
g2 = 0.1. As expected from Eq. l[65j, KSB, KSBtr and KSBl have 
the same limit for large PSF (A — >■ 1), while the bias for KSB3 
is the largest in the limit of a very wide PSF. As noted before, the 
PSF correction in all KSB variants introduces a negative bias which 
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Space based 
Ground based 



Figure 6. Dependence of A as defined in Eq. (59) on the observed elHptic- 
ity. Red line represents the case of PSF with FWHM = 0.5Re and /3 = 2 
to mimic a space-based observation and the green line the case of a PSF 
vi'ith FWHM = 5Re and /3 = 5 to mimic a ground-based observation. In 
the case of a infintely wide PSF we would have A/A(0) = 1 independently 
on the ellipticity. 



partly compensates (or even overcompensates) the overestimate by 
KSB and KSBtr from the lens mapping. Since KSB3 is essentially 
unbiased for unconvolved ellipticities, any PSF correction neces- 
sarily lowers the shear estimate. 

Finally, we introduce the weight function into the moment 
measurement and study the response of the four methods in this 
situation. The result is shown in Fig. [TO] for a space-based (left 
panel) and a ground-based observation (right panel). For narrow 
PSFs, the methods react on weighting as in the previous Section, 
where the PSF was neglected (see Fig. while the response is 
milder for a wider PSF. For a narrow PSF, KSBtr and KSB3 are 
essentially unbiased, and KSBtr remains fairly unbiased when the 
PSF width increases. From the comparison between Figs.|9land|10| 
we can infer the effect of weighting on the shear estimates. The bi- 
ases of most methods are lowered because the ellipticity of the con- 
volved source is lower, hence a circular weight function does not 
significantly affect the ellipticity measurement. However, in partic- 
ular KSBl shows concerning dependence on both the presence of 
a weighting function and the width of the PSF: Even though KSBl 
seems fairly unbiased in the right panel of Fig. (TO) other values of 
the width of the weight function would lead to less optimal results. 

We are aware that our tests are of somewhat approximate na- 
ture in the sense that the characteristics of the simulated images 
only coarsely resemble that of realistic survey data. The real-life 
performance of all KSB variants will depend on peculiar proper- 
ties of the surveys to be analyzed, such as the shape of the PSF, the 
depth of the observation, etc. However, two findings from our result 
can be considered robust: KSB3 shows the least amount of bias and 
the weakest dependence on the width of the weighting function, as 
long as the PSF remains narrow with respect to the galaxy size. KS- 
Btr has a more pronounced dependence on the weighting function, 
but reacts only weakly on changes of the PSF width. 



5 CONCLUSIONS 

We have assessed the assumptions underlying the KSB method 
for measuring gravitational shear from the images of ensembles 




Figure 7. Dependence of A as defined in Eq. (59) on the size of the PSF for 
a fixed value of the observed ellipticity. 




KSB 
KSBtr 
KSB3 
KSBl 



2 4 6 8 10 12 14 16 18 2C 
PSF size [Re] 

Figure 8. Shear estimate gi as a function of the PSF size for a Sersic-type 
galaxy image as provided by KSB (red line), KSBtr (green line), KSB3 
(blue line), and KSB 1 (magenta line) for a fixed value of the preconvolved 
ellipticity corresponding to g = (0.4, 0.1)). 



of lensed galaxies. KSB has the great advantage of being model- 
independent since it expresses the lensing-induced shape change 
by a combination of moments of the surface-brightness distribu- 
tion. However, several assumptions underlying the derivation of the 
method and its practical implementations turn out to be violated 
more or less severely in realistic situations. We can summarise our 
results as follows: 

(i) KSB defines a shear estimate for each individual galaxy, de- 
fined as the shear that would describe the observed ellipticity if 
the object was perfectly circular prior to lensing. In other words, 
it is assumed that the instrinsic ellipticity of the individual object 
vanishes. The true shear is then computed averaging these shear 
estimates within a region where g is assumed to be constant. This 
is in general not equivalent to averaging the ellipticities of each 
individual object and then computing the true shear: averaging ob- 
served galaxy ellipticities and measuring the shear do not commute 
because not the individual intrinsic ellipticities can be assumed to 
vanish, but only their average. We show that the difference between 
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Figure 9. Shear estimate g\ as a function of the applied shear for noise-free Sersic-type galaxy images as provided by KSB (red line), KSBtr (green 
line), KSB3 (blue line), and KSBl (magenta hne). In the left panel we choose a PSF with FWHM = 0.5-Re and /3 = 2 to mimic a space-based 
observation, while in the right panel we choose FWHM = 5/?^ and /3 = 5 to mimic a ground-based observation. No weight function has been used 
to compute moments. 
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Figure 10. As in Fig.|9]but employing a weighting function in the moments computation. The width of the weighting function was set to the apparent 
size of the objects, a = \/Tr(Q) 



the two approaches is a function of the variance of the intrinsic el- 
Upticity distribution. The error introduced this way depends on the 
variants of KSB used, the size of the PSF, and the width of the 
weight function. It is normally in the percent range. 

(ii) The definition of the KSB shear estimate relies on the as- 
sumption that the shear is small. However, this is only true after 
averaging. For a single object, the reduced-shear estimate g is of 
the same order as the ellipticity x- This leads to a relation between 
g and x which is correct only to first order. This situation can be 
improved considering only linear terms in X9 in ths derivation of 
pah (j^sgi)^ Qj. considering consistently terms up to third order in 
Xg (KSB3). We also show that the approximation of P"'' by half 
of its trace (KSBtr), although not mathematically justified, yields a 
better g-x relation compared to KSB. 

(iii) KSB, KSB 1 and KSBtr in absence of PSF convolution tend 
to overestimate the shear, while KSB3 gives an almost perfect re- 
sult. 

(iv) KSB and KSB 1 depend strongly on the width of the weight 



function used in the moment measurements, while KSBtr is more 
robust and KSB3 is almost independent of it. 

(v) KSB does not perform any PSF deconvolution, but gives 
only an approximate correction for the effects of the PSF. We show 
that this correction would be equivalent to a proper deconvolution 
from a circular PSF only in the case of a circular source, otherwise 
the improper PSF correction lowers the shear estimate. 

(vi) The overestimate due to the wrong relation between g and 
X and the underestimate due to the inappropriate PSF correction 
tend to compensate each other. For a narrow PSF (space-based ob- 
servation), KSB3 is the variant with the least bias, while KSBtr is 
the best method for wider PSFs (ground-based observation). 

(vii) The choice of the width a of the weight function could 
be utilized to reduce the measurement bias. In principle, a can be 
tuned according to the size of the PSF and to the galaxy ellipicities 
such that the shear estimate ends up to be almost unbiased. How- 
ever, practically this is only feasible for the galaxy ensemble as a 
whole, whereas choosing a such that shear estimates are unbiased 
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for each individual galaxy is of similar difficulty as estimating the 
shear. 

(viii) KSB can correct only small anisotropics in the PSF {q <^ 
1). It is not possible to extend the formalism to allow more precise 
corrections since that would imply the calculation of derivatives of 
the observed surface brightness, which is not feasible since 7°''* 
is a noisy quantity. An improper correction of the PSF anisotropy 
introduces a bias which propagates to the final measurement of the 
shear in an almost unpredictable way. 



Kii — (Qiiiiii — <3iiii22 ~ Q112222 + Q222222) (A-17) 

K12 = K21 = 2((5"iiii2 — (5122222) (A18) 

K22 = 4(Q'i'ni22 + Q'1'12222) (A19) 

-Dill = Qiiiiii — 3Q111122 + 3(5iiii22 — Q222222 (A20) 



APPENDIX A: 

In this Appendix, we list expressions for the tensors defined in the 
paper in terms of moments of the light distribution: 



- 2xi + 2 ^ + 2 



Tr{Q) Tr{Q) 



P22 = -2#^ - 2X2 + 8;^% + 2 



Tr{Q) 



Tr{Q) 



Tr{Q) Tr{Q) Tr{Q) Tr{Q) 



Tr{Q) Tr{Q) Tr{Q) 



Tr{Q) Tr{Q) Tr{Q) Tr{Q) 



Tr{Q) Tr{Q) Tr{Q) 



Tr(Q) Tr{Q) Tr{Q) 



R,,,=2^^+4^^-2,^''' 



Tr{Q) Tr{Q) Tr{Q) 



Tr{Q) Tr{Q) Tr{Q) Tr(Q) 



R222 = 2 



K22X2 , ,c'5'll22X2 ic*?!!!: 

— + lo TTT^ Id 

tr{Q) 



+ 16 



tr{Q) Tr{Q) Tr{Q) 



Li = Qiiu — Q2222 
1/2 = 2(Q'iii2 + (52221) 

Bll = Qllll — 2(5ll22 + (52222 

S12 = S21 = 2((5'ni2 - Q'1222) 
-B22 = 4(5n22 



(Al) 

(A2) 

(A3) 

(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 

(A12) 
(A13) 
(A14) 
(A15) 
(A16) 



i?112 = D121 = -D21I = 

= 2((5lllll2 — 2(5lll222 + (5222221) 

D122 = D212 = D22I = 4((5llll22 ~ <5ll2222) 

-D222 = 8(5'l 11222 

Chilli = (5llllllll ~ 4(5iiiiii22 + 6(5llll2222 

- 4Q'i' 

1222222 + (522222222 

t/2111 = 2((5lllllll2 — 3(5lllll222 
+ 3(5lll22222 ~ (512222222) 

[^1211 = J/ll21 = J/lll2 = t/2111 
[^2211 = 4((5llllll22 ~ 2(5llll2222 + (5ll222222) 



(A21) 
(A22) 
(A23) 

(A24) 

(A25) 
(A26) 
(A27) 



C/2112 = C/2121 = C/1122 = C/1221 = C/1212 = C/2211 (A28) 

C/2221 = 8((5'l"lll222 — <5'l"l22222) (A29) 

[^2122 = C/2212 = t/1222 = 1/2221 (A30) 

t/2222 = 16(5nii2222 (A31) 

Note that the moments Q'ij,,k have to be computed using the first 
derivative of the weight function with respect to 6^. Similarly, 
Qij,,k must be computed using the second derivatives of the weight 
function. 
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